Orbital contribution to the magnetic properties of nanowires: Is the orbital 

polarization ansatz justified? 
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We show that considerable orbital magnetic moments and magneto-crystalline anisotropy ener- 
gies are obtained for a Fe monatomic wire described in a tight-binding method with intra-atomic 
electronic interactions treated in a full Hartree Fock (HF) decoupling scheme. Even-though the use 
of the orbital polarization ansatz with simplified Hamiltonians leads to fairly good results when the 
spin magnetization is saturated this is not the case of unsaturated systems. We conclude that the 
full HF scheme is necessary to investigate low dimensional systems. 
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In the bulk of ferromagnetic transition metals it is 
well known that the orbital magnetic moment L is 
quenched and that the magneto-crystalline anisotropy 
energy (MAE) is very small as a result of strong electron 
derealization and crystal field effects compared with 
those of intra-atomic Coulomb interactions. In nano- 
objects the dimensionality or coordination is reduced so 
that the influence of these interactions, responsible for 
Hund's rules in the free atom, becomes more and more 
important and both the spin and orbital magnetic mo- 
ments increase dramatically. This is seen in experiments 
on chains of Co atoms at step edges of Pt(997)|lj and 
Co single atoms or nanoparticles deposited on Pt(lll) 
in which orbital moments as large as l.I/is per atom 
have been measured 2] , associated with a considerable 
enhancement of the MAE. 

On the theoretical side, in the Local Spin Density Ap- 
proximation (LSDA) or in simplified tight-binding (TB) 
Hartree-Fock (HF) schemes the intra-atomic Coulomb in- 
teractions are treated in an average manner so that the 
distribution of electrons between the orbital states of op- 
posite magnetic quantum numbers m is poorly described, 
especially in low dimensional systems. As a result these 
approximations yield underestimated values of L, even 
though these values increase when the dimensionality is 
lowered. Eriksson et ai.|3| have proposed to correct for 
this effect by adding a term proportional to — L 2 /2 in 
the Hamiltonian, treated in mean-field, which will be re- 
ferred to as Orbital Polarization Ansatz (OPA) in the 
following. The effect of this term is obviously to in- 
crease < L >. A more rigorous way of obtaining both 
the spin and orbital moments is to solve the HF equa- 
tions by taking into account all intra-atomic terms in the 
decoupling with all matrix elements of the Coulomb in- 
teraction C/ 7l727374 = (7i(^),72(r > ')|pz^|l73(T , ),74(r / )}, 
where 7$ are atomic orbitals, expressed in terms of the 
three Racah parameters A, B and C, for d electrons 
and a system of homonuclear atoms. Starting from this 
Hamiltonian Solovyev et al. have shown, in an elegant 
work, that the OPA cannot be derived analytically from 



the HF Hamiltonian except in some very special cases 
and that, even in the latter, the proportionality factor 
is not B as usually assumed but 3B/2. Very recently 
Nicolas et al. have discussed the effect of orbital polar- 
ization, using either a Stoner-like TB Hamiltonian with 
the OPA or an HF Hamiltonian in which the one and 
two orbital matrix elements of the Coulomb interaction 
are treated exactly in the spherical harmonics (SH) basis 
but three and four orbital terms are neglected. These 
latter terms depend both on B and C in the SH basis 
which results in a symmetry breaking that they claim to 
overcome by averaging over different orbital basis. On 
the opposite, a recent work by Xiangang Wan et ai.Q is 
based on a complete HF decoupling. However their ef- 
fective intra-atomic potential (see Eq.4 of their work) is 
the same as in LSDA+U while the TB part of their total 
Hamiltonian is not spin polarized. As a result when the 
approximations leading to the Stoner model are carried 
out in their Eq.4, it does not lead to the correct Stoner 
parameter. 

It is thus of fundamental importance to investigate 
the ability of the full HF scheme to predict large L and 
MAE in nano-objects, and to check whether the OPA 
can account for these effects. In this paper we com- 
pare, on the simple model of a monatomic wire, the 
results given by the full HF decoupling and two cur- 
rently used simplified Hamiltonians corrected or not by 
the OPA term. We use a TB model in a minimal orthog- 
onal basis set of d valence orbitals \i, 7, a) = <8> |cr), 
of spin a and orbital 7 centered at site i. In the fol- 
lowing 7 will either denote cubic harmonics (CH) (7 = 
A = d xy , d yz , d zx , d x 2_ y 2, d 3z 2_ r 2) or spherical harmonics 
(7 = m = —2, —1, 0, 1, 2). Our Hamiltonian H can be ex- 
pressed as the sum of a standard one-body TB Hamilto- 
nian Hq (determined by the bare d level eqj an d hopping 
integrals) and an electron-electron interaction Hamilto- 
nian Hi nt in which only on-site electron-electron interac- 
tions are considered. The standard Hartree Fock decou- 
pling leads to the one-electron Hamiltonian (denoted as 
HF1) which, in the second quantization formalism can be 
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written 

•^int = (^74727371 ( C l'y4a C n3c) c i-/ 2 a' C Hl^' 

171727374 
crcr' 

— C/74727173 ( C l 74 o- C i730-') C i 72( T' Ci 7lO'^ (1) 

The expressions of the matrix elements l^ 7i727 3 7 4 obvi- 
ously depend on the atomic basis, but the resolution of 
the full Hartree-Fock Hamiltonian (namely without any 
approximation) must lead to the same results whatever 
the basis. However, the use of CH is quite attractive for 
discussing the OPA since in this basis the three and four 
orbital matrix elements of the electron-electron interac- 
tion are proportional to the Racah parameter B only^jj. 
Moreover in CH the different values of the two orbital 
matrix elements Ux^x^ and Ux^^x (A ^ /i) only dif- 
fer by terms proportional to B. The average values: 

(1/4) E MlW 4A ^AmAm and (1/4) U Xm* are im- 

pendent of A and are given by U = A — B + C and 
J = 5-E>/2 + C§| while the one orbital terms Uxxxx are 
all equal to U + 2 J. This leads us to define (U, J, B) as 
a new set of parameters. The two orbital terms Ux^Xij. 
(resp. Uxfifix) can then be expressed in terms of U and 
B (resp. J and B) while the three and four orbital terms 
are proportional to B only. As already stated, this is no 
longer true in the SH basis. 

When B is neglected in the above Hamiltonian HF1, 
we recover the model (hereafter referred to as HF2) that 
has been used in our previous studies 9] (Ux^x^ — U and 
UxnnX — J f° r an y Pah of different orbitals A and /1 
and no three and four orbital terms) in which spin-flip 
terms were omitted since the spin-orbit coupling inter- 
action was not taken into account. Starting from this 
Hamiltonian, keeping only the diagonal terms and replac- 
ing each orbital population of a given spin by its average 
value, leads to a Stoner-like Hamiltonian (called HF3) 
that we have also investigated since it has widely been 
used in the literature |Toj: 

^f 3 = YJVrtNi - alMiWcl^dxo . (2) 

iX,a 

In this hamiltonian I — (U + 6 J)/5 is the Stoner param- 
eter while Ni and Mj are, respectively, the total charge 
and moment on site i. U e g is equal to (%U — 2J)/10 if one 
derives HF3 from HF2 as explained above. Since here we 
are interested in systems with geometrically equivalent 
atoms (i.e., 2Vj = N,Mi — M) we can choose the energy 
zero in all hamiltonians as eo+U e &N so that the first term 
in Eq.(3) disappears from the total hamiltonian HF3. 

The spin magnetism is governed by the Stoner param- 
eter I that will be kept constant in all our calculations 
and determined so that it reproduces the experimental 
value in the bulk bec phase. 



From the above discussion it is clear that HF2 differs 
from HF1 by terms proportional to B, this is also true for 
HF3 as far as this hamiltonian is justified. Eriksson et 
alQ have proposed to introduce an OPA term to account 
for this difference. This term is written in meanficld 
AE OP = -f^Ez^) 2 which reduces to -^^(L^) 2 
when the spin and orbital moments have the same quanti- 
zation axis z (which is strictly verified along high symme- 
try directions). The corresponding Hamiltonian is then: 

H OP = -B{L Z ) ^[L^yct^c^v (3) 

where [L z ] 7 y are the matrix elements of the local or- 
bital moment operator Li Z , [L z ] 77 / is spin independent 
and diagonal in the SH basis when the orbital momen- 
tum quantization axis of the SH orbitals is rotated so 
that it coincides with the spin quantization axis. This is 
no longer true if the SH orbital momentum axis is along 
a crystallographic axis which is not parallel to the spin 
quantization axis, or when [Z z ] 7 y is expressed in the CH 
basis. Finally the last term of our Hamiltonian takes 
into account the intra-atomic spin-orbit interactions de- 
termined by the spin-orbit coupling parameter £. 

A monatomic wire of a transition metal is a handy 
system to compare the results given by the various mod- 
els described above. The parameters of the model are 
chosen to mimic Fe which is assumed to have N — 7 va- 
lence d electrons per atom in the bulk as well as in the 
wire. The hopping integrals dda, ddir and ddS are cho- 
sen proportional to (-6,4,-1) and decrease with the inter- 
atomic distance according to a R~ 5 law. The numerical 
value of dda is fitted to the bulk d band width of Fe 
(Wd = 6eV) which leads to dda = -0.749eV at the bulk 
nearest neighbor distance (d=4.7a.u.). First and second 
nearest neighbor hopping integrals have been taken into 
account. The Stoner parameter is / = 0.67eV. The spin- 
orbit coupling parameter is taken from a previous work 
(£ = 0.06eV)|ll|. It is well known that the parame- 
ter U is strongly screened in metals. In particular in a 
recent paper Solovyev^ has shown that this parame- 
ter is almost independent of the bare interaction. From 
Fi g.l of this reference it can be deduced that U ~ J in 
Fe 13]. In that case I = 7 J/5 so that U = J = 0.48eV, 
a numerical value in good agreement with that given by 
Solovyev. Finally, as in previous works^j, we have taken 
S = 0.14j|3. 

When applied to bulk Fe, the complete HF decoupling 
yields (2S Z ) = 2.12// fl and (L z ) = 0.08^ B when B = 
(HF2 model), and (2S Z ) = 2.11/x B and (L z ) = 0.12/i B 
when B is taken into account (HF1 model). Then, we 
have compared the results derived from the five models 
(HF1, HF2 and HF3 with and without H p) for the spin 
and orbital moments with magnetizations along the wire 
(8 = 0) and perpendicular to it (6 = n/2) and the corre- 
sponding magnetocrystalline anisotropy energy (MAE) 
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1.45 


0.37 
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0.37 


1.31 
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0.49 


0.25 
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0.60 
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(25 z (0)} 


1.51 


1.24 


1.23 


0.94 


0.78 


(2S z (n/2)) 


1.51 


1.23 


1.24 


0.93 


0.94 


(i.(0)> 


0.33 


0.19 


0.39 


0.24 


1.07 


(L.(tt/2)) 


0.21 


0.10 


0.18 


0.08 


0.15 


MAE 


-0.7 


-0.3 


1.5 


0.0 


6.2 



TABLE I: The spin ((2&)) and orbital ((L z )) magnetic mo- 
ments (in hb per atom) for a monatomic Fe wire and two 
magnetization orientations (parallel (9 — 0) and perpendic- 
ular (9 = 7r/2) to the wire) and the corresponding magne- 
tocrystalline anisotropy MAE (Etat(ir/2) — E to t(0)) in meV 
per atom for two interatomic distances. 



AE = E tot {6 = tt/2) - E tot {6 = 0) where E tot is the 
total energy per atom of the system. Two interatomic 
distances have been considered: the bulk interatomic dis- 
tance at which the spin magnetization is saturated and a 
shorter distance (4.25a.u.) corresponding to unsaturated 
spin moments. The results are given in Table 1. 

Let us first discuss the wire at the bulk interatomic 
distance. All models agree to predict saturated spin 
magnetization, i.e., the spin magnetic moment is 3fiB 
to less than a few 10 _3 /is. As a consequence the effec- 
tive atomic orbital levels with down spin are identical in 
HF2 and HF3 models since U = J. This is no longer 
true for the up spin orbitals for which the atomic lev- 
els are orbital dependent with HF2 and not with HF3. 
However the average atomic level is the same in both 
models. Therefore the orbital moment, which arises only 
from the spin down band, the spin up band being filled, 
is almost identical in both models similarly to the total 
energy (see Table 1). As expected the orbital moments 
for both magnetization orientations and the associated 
MAE, even though reinforced compared to the bulk ones, 
are largely underestimated by the HF2 and HF3 models 
with B = compared to those predicted by the complete 
HF decoupling (HF1). When the OPA term is added to 
the HF2 and HF3 hamiltonians, the results given by the 
latter models become in fair agreement with those ob- 
tained from HF1 for the orbital moment while the MAE 
is well reproduced. 

The above trends completely change when the inter- 
atomic distance is shortened to 4.25a.u.. It is first seen 
that the spin moment depends on the model. In this re- 
spect the HF2 model is much better than the HF3 one. 



Moreover, taking into account the OPA term leads to an 
increase of the orbital moments for both magnetization 
orientations which are rather close to the HF1 results for 
HF2 but not for HF3. 
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FIG. 1: (L z ) and MAE as a function of B/J from HF1 and 
HF2 (HF3 results are undistinguishable from the HF2 ones) 
for a magnetic Fe monatomic wire (d =4.7a.u.). 

To summarize this discussion we can state that the 
OPA is rather good for saturated spin magnetization 
while for the unsaturated case it leads to results depend- 
ing critically on the approximations made concerning the 
electron-electron interaction hamiltonian. In order to 
verify that the good performance of the OPA for the 
saturated spin magnetization is not due to the partic- 
ular value of B, we have studied the variation of (L z ) 
at 9 — and 9 — tt/2 and the associated MAE as a 
function of the ratio B/J. The results (Fig^l show that 
the OPA gives the right trends on the full domain of 
B/J values that we have investigated. In particular an 
abrupt variation of (L z ) at 9 = occurs around a critical 
value of B/J ~ 0.09 above which the upper 8 band (the 
corresponding eigenfunctions being mostly linear combi- 
nations of SH with \m\ — 2) of minority spin becomes 
empty. 

Even if the OPA works reasonably in the saturated spin 
magnetization case for determining (L z ) and the MAE, 
this does not mean that it reproduces the band structure 
correctly. Let us first note that for 9 = 0, the eigen- 
functions have a largely dominating single SH character 
while at 9 = tt/2 they are almost pure single CH orbitals. 
The band structures corresponding to HF1 and HF2 are 
drawn in Fig (the band structure of HF3 is close to that 
of HF2). At first sight they look quite similar. However a 
closer examination reveals some differences. Let us first 
comment on the majority spin bands at 9 = 0. While 
the splittings of the |m| = 2 (S) and \m\ = 1 (tt) bands 
are respectively given by 2£ and £ with the HF1 model, 
they become 2£ - 4B{L Z ) and £ - 2B(L Z ) with both the 
HF2 and HF3 models, respectively. In addition the m 
character of the bands is reversed, i.e., the m=2(l) band 
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FIG. 2: HF1 (top) and HF2 (bottom) band structure (re- 
ferred to the Fermi level) for a magnetic Fe monatomic wire 
{d =4.7au) with a magnetization parallel (9 = 0) and per- 
pendicular (6 — tt/2) to the wire. All results are obtained 
for U/J — 1 and B — 0.14 J save for the dotted band struc- 
ture of the the top right panel obtained for U/J — 1.34 and 
B = 0.14J. 

is above the m=-2(-l) band in the HF1 while it is the 
opposite with the HF2 and HF3 models. This inversion 
does not occur in the minority spin bands and the split- 
tings of the 5 and n bands are not exactly the same with 
the HF2 and HF3 as with the HF1 models. At 9 = tt/2 
all models agree that for U = J there are almost no band 
splittings and that the removals of degeneracy around the 
midpoint between T and X are more pronounced in the 
minority bands than in the majority ones. 

Finally it is interesting to study the variation of (L z ) 
and of the MAE with the HF1 model when the ratio U/ J 
is varied by keeping the Stoner parameter fixed. Indeed 
this ratio is not perfectly known. The results are shown in 
FigEH Abrupt variations of (L z ) are observed at U/J~ 
1.34 when 9 = tt/2 and U/J ~ 3.25 when 9 = 0. They 
correspond respectively to the occurrence of a splitting of 
the 5 bands (see Fig|2J) and to the complete filling of the 
lowest 5 band of minority spin. These abrupt changes of 
(L z ) are associated with a change of sign of the variation 
of the MAE as a function of U/J. 

In conclusion we have studied orbital polarization ef- 
fects for a Fe monatomic wire with various HF Hamilto- 
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FIG. 3: (L z > and MAE as a function of U/J from HF1 for a 
magnetic Fe monatomic wire (d =4.7au). 



nians in a tight-binding scheme: a full HF Hamiltonian 
(HF1) including all the Coulomb interaction matrix ele- 
ments, a simplified one (HF2) neglecting the Racah pa- 
rameter B, and finally a Stoner-likc Hamiltonian (HF3). 
OPA has then been reintroduced in HF2 and HF3 as 
proposed by Eriksson et ajQ. With HF1 we predict that 
very large values of L and MAE are possible in agree- 
ment with existing experiments. The same trends are 
obtained by adding the OPA to simplified Hamiltonians 
when the spin moment is saturated, however noticeable 
differences appear in the band structure since some split- 
ting and band characters are wrongly reproduced. This 
fair agreement strongly deteriorates when dealing with 
an unsaturated system, especially with the Stoner-like 
model. It is thus of prime importance to use the HF1 
model for the study of systems with much more com- 
plex geometries (surfaces, clusters, break junctions), in a 
realistic s, p and d basis set, or to implement it in ab- 
initio codes. Indeed from our results giant anisotropy of 
magneto-resistance in low dimensional systems such as 
break junctions is expected^). 
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